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Au+Au, s^/^ = 200 A GeV measurements at RHIC, obtained with the PHENIX, STAR, PHOBOS 
and BRAHMS detectors, have all indicated a suppression of production, relative to an appropriately 
normalized NN level. For central collisions and vanishing pseudo-rapidity these experiments exhibit 
suppression in charged meson production, especially at medium to large transverse momenta. In 
the PHENIX experiment similar behavior has been reported for tt" spectra. 

In a recent work [J] on the simpler D+Au interaction, to be considered perhaps as a tune-up for 
Au-|-Au, we reported on a pre-hadronic cascade mechanism which can explain the mixed observation 
of moderately reduced p± suppression at higher pseudo-rapidity as well as the Cronin enhancement 
at mid-rapidity. Here we present the extension of this work to the more massive ion-ion collisions. 

Our major thesis is that much of the suppression is generated in a late stage cascade of colourless 
pre-hadrons produced after an initial short-lived coloured phase. We present a pQCD argument 
to justify this approach and to estimate the time duration Tp of this initial phase. Of essential 
importance is the brevity in time of the coloured phase existence relative to that of the strongly 
interacting pre-hadron phase, the latter essentially an interactive cascade. These distinctions in 
phase are of course not strict, but adequate for treating the suppression of moderate and high p± 
mesons. 

PACS numbers: 25.75.-q,12.38.Mh,24.10. Jv 



I. INTRODUCTION 



The specific question we deal with in this work is the suppression of medium to high transverse momentum meson 
yields observed in the experimental measurements 0, d, for Au-|-Au at 130 GeV and 200 GeV. The experiments 
have focused on the rj and pj^-dependence of the ratio of these yields to those in p+p, for charged particles: 



R AA/NN = ^— J, '\ 1 

where ^coii is a calculated number of binary NN collisions occurring in Au+Au at a designated energy and centrality. 
One can also, of course, just compare the yields directly to the data without reference to ratios. 

The simulation code LUCIFER, developed for high energy heavy-ion collisions, has previously been applied to both 
SPS energies s^/^ =(17.2,20) A GeV [1,0, H and to RHIC energies s^/^ = (56,130,200) A GeV [1, l|. Although 
nominally intended for dealing with soft, low p_L interactions, it is possible to introduce high p±^ hadron spectra via 
the NN inputs, which form the building blocks of the simulations, and thus to examine the effect of re-scattering, and 
concomitant energy loss, on the high spectra 0]. The simulation is divided into two phases I and II. The first 
stage sets up the participants, both mesonic and baryonic, their four momenta and positions for the commencement 
of a 'hadronic' cascade in phase II. Energy loss and interactions within what amounts to a 'hadronic' fluid in phase 
II play a key role in the eventual suppression of the transverse momentum distributions. 

This approach has already yielded a consistent description of at least the final stage phenomena leading to mea- 
surements in both D + Au [l| and Au + Au at a variety of kinematical conditions, energies and transverse momenta. 
All these results follow given only the presence at early times during the ion-ion collision of a fluid of 'pre-hadrons' 
having rather generic properties. 

Such a mechanism for the suppression of high transverse momentum jets may seem to run counter to the conventional 
pQCD view. There does exist, however, notable justification for a picture in which colourless pre-hadrons, may be 
produced at rather early times in a AA collision, and then play a key role in the further development of the system. 
First there is the work of Shuryak and Zahed |9|] on the persistence of hadron-like states above the critical temperature 
in a dense quark-gluon medium as well as similar results from lattice-driven studies on the persistence of the J/^ and 
other special hadronic states [T^. Recently a much more direct treatment has been given, on which we now mainly 
rely, by Kopcliovich, Nemchik and Schmidt |Jd|] as well as Berger [l^], who directly consider the temporal dynamics 
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of pre-hadron production from a pQCD point of view. In particular Kopeliovich et al., outline the fate of leading 
hadrons from jets, in deep inelastic scattering (DIS) on massive ions and discuss the relevance to relativistic heavy 
ion collisions. These works have the distinct advantage of relying directly on pQCD. We attempt to reproduce their 
development, as it is germane to our treatment of NN, NA and AA interactions. 

Figure (1), essentially borrowed from Refercnce(ll), displays the production of a jet, here at least a moderately 
high energy quark, in a pp or pA event and the time scales relevant to the process. Kopeliovich et al. begin by 
explicitly considering eP or eA, in which case the initiating particle in the jet creation is an off-shell photon (7*). In 
pp or AA one could equally well employ a gluon at the first vertex. The kinematics [ll| remains essentially the same 
then for nuclear induced events and we elaborate this case somewhat. 

The essential feature evident in the figure is the production of a colour neutral pre-hadron after a time tp, from 
an initially pcrturbative high pj_ quark, by what one could label as coalescence of the leading quark with an ambient 
anti-quark. This coalescence incidentally is most probable for co-moving anti-quarks, the scarcity of co-movers with 
increasing p±^ explaining the rapid fall off of even the pp meson spectrum with increasing p±. Also prominent in 
this diagram are the early perturbative gluon radiation. These radiated gluons may in due course create further pre- 
hadrons, and of course there is then a concomitant energy loss from the initial quark. Further, for leading partons, 
one expects that the hadronisation time th is considerably longer than tp, involving soft processes that eventually put 
the pre-hadron on its final hadronic mass shell. Octet gg's might also arise but will not persist or eventually coalesce 
given the repulsive or non-confining forces in such systems. 

One can summarise the kinematic arguments [ll[ most easily in the rest frame of the struck proton or nucleus. We 
use these arguments to estimate the production time tp for the pre-hadron, which will be a considerably shorter time 
Tp in the colliding frame of an AA system. For DIS on a nucleon, on its own or within a nucleus A, the production 
time in Figure (1) is estimated to be 



where Zh is the fraction of quark energy imparted to the hadron. Integrating the gluon radiation spectrum one obtains 
for the energy loss per unit length z, a time independent rate [Ulll^ rising quadratically with the hard scale Q, 



^ ^ (3) 

Uz Stt 



The colour neutralisation time is then given by: 



*p~q|(1-^'0- (4) 
The hadronisation time is related to the QCD scale factor and is usually estimated as: 

th^ .f^ (1 - Zh), (5) 

and is much longer than the colour neutralisation time, given that the QCD scale is close to 200 GeV. 

The authors of Reference (11) argue the brevity of the de-colourisation time scale is a consequence of energy 
conservation, coherence and Sudakov suppression. 

Importantly, a much reduced scale for pre-hadron creation is likely to remain true even for non-leading partons 
provided Equation (4) remains valid, i. e. when Q is comparable to Eg. Thus many of the radiated hard or even 
moderately hard gluons, in the early stages of a pp or AA collision will initiate similar pre-hadrons, and in a nuclear 
medium such large sized objects will be sufficiently numerous and of critical importance to the dynamics. One should 
keep in mind that the overall multiplicity is not large at p_L > 1 GeV/c, where the NN spectra (see Figure (5)) has 
dropped from that at the softest transverse momenta measured by close to two orders of magnitude and for Au+Au 
(Figure 7) by somewhat more. 

One concludes that the production and hadronisation processes are to some degree separate: with generally tp less 
than tf and frequently much less. Corresponding time scales in, say, the center of momentum frame for an A-l-A 
system will of course be considerably contracted. Although the colourless pre-hadron in Figure (1) is generated by an 
initially perturbative process with an anti-quark, it is the subsequent interactions with other such pre-hadrons that 
leads to the observed suppression for mesons of appreciable and even moderate transverse momentum. The pre-hadron 
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perturbatively begins life with the qq at small relative distance, and thus has a small initial size, but in light of the 
ambient momenta for such high or even moderately high energy partons, the transverse diameter rapidly increases to 
the scale of a typical hadron. Reference (11) suggests the entire growth to pre-hadronic size occurs quickly. 

One can pursue the evolution of the system of pre-hadrons via a Glauber-theory [H H based treatment of their 
interactions, or, as we do, via a standard cascade model. The end result is little different: what has been labelled 
jet suppression results. In Glauber theory the hadron-sized cross-sections produce strong absorption: and hard pre- 
hadrons simply do not remain in the final state, they are too often absorbed. In the cascade described hereafter the 
pre-hadronic medium is sharply cooled by the interactions and instead of on hard meson, many soft mesons appear, 
and at lower p±^ . 

To obtain the final cascade yields for both D-|-Au [l| and Au-|-Au, it is essential that the characteristic time tp be 
considerably less than tj. Indeed this constraint plus the appropriately large pre- hadron interaction cross-sections 
are also instrumental in generating the surprisingly large elliptical flow that has been measured at RHIG [l^ [l^ . 
References (11,12) provide a perturbative QCD justification for both such conditions, the early de-colourisation and 
the large pre-hadron cross-sections. 

For completeness we present a brief overview of the dynamics of our Monte Garlo simulation, which has already 
been described extensively in earlier works [l], 0| ■ Many other simulations of heavy ion collisions exist and these are 
frequently hybrid in nature, using say string models in the initial state [TtI. [Ts^ 19^20^21^22^2 3,_ 2^ together with 
final state hadronic collisions, while some are either purely or partly partonic 25, 26l. 27. 28l 29l IsOj in nature. Our 
approach is closest in spirit to that of RQMD and K. Gallmeister, G. Greiner, and Z. Xu [31[ as well as work by 
W. Gassing [l3,|3ll[32| and G. Wolschin [SJ]. Gertainly our results seem to parallel those of the latter authors. They 
seek to separate initial perhaps parton dominated processes from hadronic interactions occurring at some intermediate 
but not necessarily late time, or to rely on some hadronic dynamics. 

The ostensible purpose of describing such high energy collisions without explicitly introducing the parton substruc- 
ture of hadrons, at least for soft processes, was to set a baseline for judging whether deviations existed between 
simulation and experiment, which could then signal interesting phenomena. The division between soft and hard 
processes is not necessarily easy to identify in heavy ion data, although many authors believe they have accomplished 
this within a gluon-saturation picture 34, [3^ [s^. For both D+Au and Au+Au systems we separate the effects of a 
second stage cascade, a lower energy but still early hadronic cascade, from those of the first stage, a parallel rather 
than sequential treatment of initial (target)- (projectile) NN interactions. The first, quite likely partonic, stage ends 
after some average time tp. Paradoxically for us then, it is perhaps higher p±^ meson production that is most justified 
by pQGD. 

In the present work we do indeed find considerable suppression of the Au-f Au transverse momentum spectrum 
at central rapidity. The key dynamics occur presumably after the de-colourisation phase in Reference (11), i. e. 
after some time closely related to the much contracted production time in the collision center of momentum frame 

The first stage of our simulation, involving the collective interaction of the initially present nucleons, produces a 
'hot' gas or liquid of pre-hadrons which is considerably cooled in an inevitably expanding and interacting final state 
cascade. This cooling produces the observed 'jet' suppression, and not surprisingly, that suppression is appreciably 
greater for Au+Au than for D+Au [U, @]- It is clear that the early pre-hadron 'medium' resembles a rather dense 
fluid. 



II. THE SIMULATION 



A. Phase I 



Phase I of LUGIFER considers the initial interactions between the separate nucleons in the colliding ions A+B, 
but is not a cascade. It amounts to a creation of the initial conditions for the final lower energy cascade. The totality 
of events in phase I involving each projectile particle happen essentially simultaneously or one might say in parallel. 
Energy loss due to pre-hadron production as well as the creation of transverse momentum (p^), initial conditions 
for what follows, are occur in phase I. A model of NN collisions @, 0], incorporating most known inclusive cross- 
section and multiplicity data, guides phase I. The two body model, clearly an input to our simulation, is fitted to 
the elastic, single-diffractive (SD) and non-single-diffractive (NSD) aspects [s^l of high energy pp collisions [12, ISsj 
and PP data 44]. The soft QGD aspects of early A-l-B events are then modeled through their incorporation in the 
phenomenological nucleon-nucleon inputs. 

One notes it is precisely the energy dependence of the cross-sections and multiplicities of the NN input that led 
to our successful prediction ^ of the rather small (~ 13%) increase in dN'^^ /drj that was first observed between 
= and 5^/2200 A GeV in the PHOBOS data 



4 



A history of the coUisions that occur between nucleons as they move along straight hnes in phase I is recorded and 
later used to guide the determination of multiplicity. Collision driven random walk fixes the p± to be ascribed to the 
baryons at the start of phase II. The overall multiplicity, however, is subject to a modification, based as we believe 
on natural physical requirements 0]. 

The collective/parallel method of treating many NN collisions between the target and projectile is achieved by 
defining a group structure for interacting baryons, which preserves the geometry of the nucleus-nucleus collision. This 
is best illustrated by considering a prototypical proton-nucleus (pA) collision. A group is defined by spatial contiguity. 
A nucleon at some impact parameter b{x±) is imagined to collide with a corresponding 'row' of nucleons sufficiently 
close in the transverse direction to the straight line path of the proton, i. e. within a distance corresponding to the 
NN cross-section. In a nucleus-nucleus (AB) collision this procedure is generalized by making two passes: on the 
first pass one includes all nucleons from the target which come within the given transverse distance of some initial 
projectile nucleon, then on the second pass one includes for each target nucleon so chosen, all of those nucleons from 
the projectile approaching it within the same transverse distance. This totality of mutually colliding nucleons, at 
similar transverse displacements, constitute a group. The procedure partitions target and projectile nucleons into a 
set of disjoint interacting groups as well as a set of non-interacting spectators in a manner depending on the overall 
geometry of the AB collision. Clearly the largest groups in pA will, in this way, be formed for small impact parameters 
6; while for the most peripheral collisions the groups will almost always consist of only one colliding pN pair. Similar 
conclusions hold in the case of AB collisions. 

In phase II of the cascade we treat the re-scattering of the pre-hadrons produced in phase I. These pre-hadrons, 
both baryonic and mesonic in type, are not physical hadron resonances or stable particles such as appear in the 
particle data tables, which can only materialize after the hadronisation time is passed. Importantly, pre-hadrons, 
are allowed to interact starting at earlier times, viz. Tp is comparable in fact to the nominal target-projectile crossing 
time Tab ^ '2Rab/j- In practice we vary the effective time delay for pre-hadron collisions in phase II, between 0.2 
and 0.4 fm/c in the center of velocity frame. 

The pre-hadrons, when mesonic, may consist of a spatially close, loosely correlated quark and anti-quark pair, and 
are given a mass spectrum between and 1.2 GcV, with appropriately higher upper and lower mass ranges for 
pre-hadrons including strange or heavier quarks. The Monte-Carlo selection of masses is then governed by a Gaussian 
distribution. 



with mo a selected center for the pre-hadron mass distribution and w — mo/4 the width. For non-strange mesonic 
pre-hadrons mo ~ 800 MeV, and for strange too ^ 950 MeV. Small changes in too and w have little effect since the 
code is constrained to fit hadron- hadron data using Equation (6). 

Too high an upper limit for too would destroy the soft nature expected for most pre-hadron interactions when they 
finally decay into 'stable' mesons. The mesonic pre-hadrons have isospin structure corresponding to p, lu, or K* , while 
the baryons, in present calculation, are constrained to range across the octet and decuplet. 



The underlying NN interaction structure involved in phase I has been introduced in a fashion dictated by standard 
proton- proton modeling (37j . The components of this modeling, elastic, single diffractive (SD) and non-single diffrac- 
tive (NSD) were introduced above. Fits are obtained to the existing two-nucleon data over a broad range of energies 
s^^"^, using the above pre-hadrons as an intermediate stage. No re-scattering, only decay of these intermediate struc- 
tures is permitted in the purely NN calculation. Cross-sections in pre-hadron collisions were assumed to be the same 
size as hadron cross-sections, i.e. the same as the meson-baryon or meson-meson cross-section, at the same center of 
mass energy. Where hadron-hadron cross-sections are inadequately known we employ straightforward quark counting 
to estimate the scale. Specifically, meson-meson interactions are scaled to 4/9 of the known NN cross-sections, thus 
no new parameters are invoked. Indeed, since known data then constrains the pre-hadronic interaction, this approach 
is a parameter-free input to the AA dynamics. 

The explicit fit to these cross-sections is displayed in Figure (2). To achieve this representation we employed, for 
example, a parametrization of the data at high energies which actually works down to 5 GeV: 



P(m) = exp(— (m — m,^)'^ / w'^) , 



(6) 



B. Elementary Hadron-Hadron Model 



atot = 35.43 - 



33.4 

gO.545 



-t- 



45.23 

gO.458 



-f 0.308 (-3.364 -fln(s))^ 



(7) 



taken from V. Ezhela et al. [38] . 
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We have also created an extended fit to the elastic data using higher energy CDF measurements at 546 GeV and 
1800 GeV tSa]: 

= 7-0 - JSl + JSi + 0-09 (-3-364 + ln(.))^ (8) 

The second important aspect of the NN interactions are the multiplicities of produced particles, which are again 
taken from measured and analyzed data. We employ KNO [40| scaling, with very little loss of accuracy for the 
transverse momenta of leading mesons, presumably from jets; i. e. the probability of achieving a given multiplicity 
at a certain energy is simply a function of the ratio of that multiplicity to the mean multiplicity. To illustrate the 
fits, Figures (3,4) show two energies, bracketing those relevant to RHIC physics. Equally good representations of 
intervening energy multiplicity data are of course obtained. As noted by many authors 41] the KNO scaling begins 
to break down for the tails of the highest energy distribution (see Fig. 3), but in a fashion which again would have 
very little quantitative effect on our presentation. 

C. High Transverse Momenta 

A question that has yet to be addressed concerns the high p± tails included in our fits to NN interactions and 
thus in our calculations. In fact it is certainly to such high p± produced particles that the arguments of Reference 
(11) best apply. We simply include high p± meson events through inclusion in the basic hadron-hadron interaction 
which is of course an input to, rather than a result of, our simulation. Thus in Figure(5) we display the NSD 
(l/27rp_L)(d27V'=''"'^f^'^/dpidr;) from UAl (H. 

One can use a single exponential together with a power-law tail in p_L, or alternatively two exponentials, to achieve 
a fit to UAl 5^/^=200 GeV data. A sampling function of the form 

/ = p^{acM-p±/w) + &/((! + (r/a)^)), (9) 

gives a satisfactory fit to the pp data in the Monte-Carlo. 

Additionally, since we, for the moment, constrain our comparisons to the production of neutral pions in Au-|-Au 
we also present, in Figure (6), the PHENIX [i^ mid-rapidity p± yield for pp together with our representation of this 
spectrum. 

The previous theoretical discussion P m strongly suggests that even moderately high p± pre-hadrons, perhaps 
as low as 1 GeV/c, should be given hadron-like cross-sections. Furthermore, the earliest meson-meson collisions in 
our second phase cascade have an initial peak rate (in CM time) at ^ .25-. 35 fm/c, with many collisions occurring at 
considerably later times. This permits even smaller pre-hadrons to appreciably increase their transverse size before 
colliding and simultaneously suggests that most collisions occur between co-moving pre-mesons. 

D. Initial Conditions for Phase II 

The final operation in phase I is to set the initial conditions for the cascade of pre-hadrons in phase II. Again we are 
guided in this by the elementary NN data mentioned above. The energy-momentum taken from the initial baryons and 
shared among the produced pre-hadrons is established and an upper limit is placed on the multiplicity of pre-hadrons. 
A final accounting of energy sharing is imposed through an overall 4- momentum conservation requirement. 

The spatial positioning of the particles at this time could be accomplished in a variety of ways. For an AB collision 
we chose to place the pre-hadrons from each group inside a cylinder with transverse dimension determined by the 
group itself, given the initial longitudinal size of the interaction region at each impact parameter. We then allow the 
cylinder to evolve freely according to the longitudinal momentum distributions, for a fixed time Tp, defined in the 
rest frame of each group. At the end of Tp the total multiplicity of the pre-hadrons is limited so that, if given normal 
hadronic sizes appropriate to meson-meson cross-sections ^ (2/3)(47r/3)(0.6)'^ fm'^, the pre-hadrons do not overlap 
within the cylinder. Such a limitation in density is consonant with the general notion that produced hadrons can only 
exist when separated from the interaction region in which they are generated (47j . One can conclude from this that 
the pre-hadron matter acts like an incompressible fluid, viz. a liquid, a state described in the earliest calculations 
with LUCIFER [1,0]. 

Up to this point longitudinal boost invariance is completely preserved since phase I is carried out using straight line 
paths. The technique of defining the evolution time in the group rest frame is essential to minimizing residual frame 
dependence which inevitably arises in any cascade, hadronic or partonic, when transverse momentum is present. This 
dependence is due to the finite size of the colliding objects, implied by their non-zero interaction cross-section. 



6 



The collision history enumerates the number of interactions suffered by each baryon and allows us to assign baryon 
transverse momenta through random walk. Pre-hadron multiplicities, subject to the density restriction described 
above, are obtained using the KNO scaling we invoke for the elementary NN interactions and the known dependence 
on flavour. Production of baryon-antibaryon pairs is allowed, guided again by constraints from NN data. Transverse 
momentum is generated for the pre-hadrons and other produced particles, paying attention to the random walk 
increase for mesons created by multiple baryon-baryon collisions. Finally, overall energy-momentum conservation is 
imposed and with it the multi-particle phase space defined. 

III. PHASE II: FINAL STATE CASCADE 

We note that phase II is a relatively lower energy cascade with the collisions between pre-hadrons, for example, 
involving energies less than s^/^ ~ 15 GeV and even lower average energy. The sharing of energy-momentum in 
the original baryon groups with the produced mesonic pre-hadrons leads to appreciable energy loss from the initial 
baryons. This loss mirrors the initial gluon radiation energy losses described in the pQCD approach and as has 
been made clear occurs mostly before the characteristic time Tp. Phase II is as stated a standard sequential cascade 
calculation in which the pre-hadrons interact and decay as do any normal hadrons either already present or produced 
during this cascade. This cascade of course imposes exact energy-momentum conservation at the level of two body 
collisions. For Au-f Au, the effect of the pre-hadron interactions is truly appreciable, greatly increasing multiplicities 
and total transverse energy, E±^ both through production and via eventual decay into the stable meson species, but 
also, as will be seen, leading to the suppression of high p±^ hadron yields. 

We are then in a position to present results for Au-|-Au collisions at 200 A-GeV. These appear in Figures (7- 
10) for various double differential transverse momenta spectra and their ratios. Most contain comparisons with 
PHENIX 0, 13 ""^ measurements. In future work we will consider also charged particle data, where produced proton 
spectra would play a larger role with increasing p± . One expects however to see similar global behaviour for charged 
mesons. 



IV. RESULTS: COMPARISON WITH n° DATA 

Figures (5) and (6) show the input, nucleon-nucleon data for elementary production of tt^'s and charged particles. 
These transverse momentum spectra at mid-rapidity have been compared to results from PHENIX ^48*1 and UAl 'is']. 
Figure (7) contains the simulated 7r° transverse momenta spectrum for Au+Au at 77 = alongside the PHENIX 
data Q. To a large extent the suppression observed experimentally is paralleled by the simulated calculations. The 
production time Tp introduced above was given two values {2Rau/ (j)) and twice this value. The variation with this 
initial state time, a parameter in our modeling, was small. (7)'s here is the longitudinal Lorcntz factor defined above 
and introduced for each baryon group separately in its rest frame. 

It's evident from the figures that, most emphatically, one cannot ignore final state cascading. Moreover, with the 
assumptions we have made, the most crucial being perhaps the early commencement time for such cascading, the 
suppression cannot be considered as necessarily a sign for production of a quark-gluon plasma: the results may perhaps 
equally well signify only the creation of a pre-hadron dominated medium after a short initial delay Tp. We repeat: 
despite the apparently short time Tp at say s^^^ =200 GeV, one finds in practice that the phase II collisions have 
an effective production or delay time of 0.25 — 0.35 fm/c and continue for some tens of fm/c. 'Co-moving' collisions 
dominate the cascade. 

It's instructive to dissect the contributions of the different phases of the simulation, i.e. to separately show the 
spectrum at the conclusion of phase I from that which results after both phases I and II are complete. In Figure (8) 
the TT^ transverse momentum yield is shown for both these cases against the experimental data. It is immediately 
evident that the many virtual NN collisions in phase I produce a much elevated p± output and that this is in turn 
reduced by more than an order of magnitude by collisions with other pre-hadrons in phase II. Part of this effect is 
through inclusion in the dynamics of at least a kinematical treatment of energy loss. Thus the above analogy of an 
initial hot gas cooled by final state interactions during expansion, is apt. 

One might well turn this around and declare that the final state scattering of a given pre-hadron with its co-movers 
has cut down the Cronin effect, a reduction which suggests the applicability of the term 'jet suppression' through 
final state interactions. One notes parenthetically that particles lost at high p± are compensated for by an increase 
in those appearing at the lowest p±^s. This is of course part of the effective cooling observed. 

A second and equally important criterion for the simulation is the maximum density imposed on initial conditions 
for phase II. Figure (9) casts some light on this and on another issue, the actual transverse energy density attributed 
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to the earliest stages of the collision. We have included in this figure the charged dN/d-q spectra, including a BRAHMS 
charged meson result (compared to a simple fit): 

(a) for the totality of 'stable' mesons in phase I+ll 

(b) the same result for phase 1 alone when only decays of pre-hadrons but no phase II re-interactions are permitted, 

(c) pre-hadrons in phase I with no decays 

It is evident that some 2/3 of the summed tranverse energy E^^ is generated in the second expanding phase II 
when the system is increasing both longitudinally and transversely. The initial, very early, generated is then 
reduced commensurately, falls well below the Bjorken limit and is hence not all available for initial coloured 'plasma' 
generation. In present calculations at the initialisation of II, and keeping in mind the average masses assigned to 
pre-hadrons centered at 0.8 to 1.0 GeV, the ambient transverse energy densities are < 1.8 GeN / fw?) for the shortest 
initial time Tp chosen and correspondingly less for longer times. Lattice calculations [l3|, however, do consider such 
densities adequate for observation of the so-called cross-over transition. 

A. Ratios of Au+Au to NN Production 

In Figure (10) we display the often discussed ratio 0, d, d, [13] introduced in Eqn.(l). The discrepancies between 
simulation and data are magnified but the general agreement shown in Figure (3) remains intact. The direct double- 
differential cross-sections fall many orders of magnitude over the measured range of but the ratios to normalized 
proton- proton data cover a much com pre ssed range of p^. In this figure we presented PHENIX 7r° data from an 
early paper but more recent data [5l| tells a somewhat modified tale. The dip in R[AuAu/NN] at the lowest 
transverse momentum point pj.= 1.25 GeV/c, implying a maximum in this ratio at slightly higher p_L, has apparently 
disappeared. This is so for both the — 5% and — 10% central data resulting in curves in better agreement with the 
corresponding theoretical calculations at the measured p/perp's. The equivalent charged particle ratios [3, H, d, [H^l 
still exhibit the apparent maximum, but in a region where the proton spectra show considerable activity. One might 
also question the use of the ratio R[AA/NN] at very low p^, involving as it does the number of binary collisions. A 
relevant or preferred divisor there might be the number of participants. 

The pQCD arguments given above for the early de-colourisation of the initial medium would only apply to sufficiently 
hard mesons, perhaps p^_ ~ 1 GeV/c or greater. Nevertheless, our soft hadronic approach for the lower transverse 
momenta does apparently give a good description of the spectra at the lowest p±Js as well. 

In Figure (10) we present LUCIFER calculations for 0% — 10% centrality, but the shape is generic at least for 
reasonably central collisions. As one approaches extreme peripherality the ratio will at all p±^ eventually approach 
unity, provided of course the number of binary collisions Ncoii is also adjusted. 

It's worth repeating that no additional adjustable parameters, beyond the times for production, Tp, and hadroni- 
sation (decay of pre-hadrons), r/, were, nor would such an approach be consistent with our goals. We seek only a 
qualitative and reasonably quantitative understanding of the observations. 

The simulations performed here overestimate, but not by very much, the suppression of the spectrum at the lowest 
Pj^, but yield an overall, perhaps surprisingly, accurate description. 

V. CONCLUSIONS 

An alternative interpretation of the high transverse momentum suppression in Au-|-Au has been presented. We 
have employed a pQCD argument to justify the production of pre-hadrons at early times within our simulation, 
and the presence of relatively large objects leads to a great deal of scattering in the final state cascade. It's this 
rescattering that produces the 'jet' suppression in our picture. Of course, and paradoxically, it is to the higher pj_ 
mesons at RHIC that the pQCD argument best applies. It is, however, interesting to extend our treatment to the 
lower reaches of the spectrum for the ratio R[AuAu/NN]. Certainly there is still a silent elephant lurking in the 
background: the observation of rather large elliptical flow in the meson spectrum [l6j . These flow measurements are 
most easily reproduced theoretically in parton cascade models when hadron-sized cross-sections are invoked jisl , and 
cross-sections of this magnitude are what our pre-hadron cascade employs. 

In the work on D-|- Au plj] the use of such a pre-hadron spectrum exposed most clearly the simple role of dynamically- 
driven geometry in ratios of BRAHMS p± spectra at varying ry, suppression in this case actually occurring only at 
high rj. 
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For Au+Au the modification of the p± spectra during phase II of LUCIFER is considerably increased and suppres- 
sion is seen even at 77 = 0. Certainly, the RHIC experiments are creating unusual nuclear matter, at high hadronic 
and energy density, with exciting consequences. Surely more experimental exploration is required. 

One important, if controversial, aspect of our combined pQCD justification of pre-hadrons together with their 
subsequent phase II cascade, is the likelihood that hadronic measurements are not necessarily probing the earliest 
coloured stages of the AA interaction. To illuminate this coloured phase it may be necessary to invoke dilepton or 
direct 7 signals. In fact 10 GeV jets are, for the full RHIC energy s^/^ = 200 GeV, not perhaps the 'purest' jets in 
a pQCD approach, and one can perhaps reach even higher transverse momentum, p± > 10 — 12 GeV/c, for signs of 
deviation from the so far observed suppression. 
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FIG. 1: Schematic drawing, borrowed from Kopeliovich et. al [Tl|], of the perturbative formation of a pre-hadron from an 
off-sheU gluon g* incident on a quark in the rest frame of nucleus B in an A+B or P+B event. The quark radiates gluons and 
eventually combines with a perturbatively produced anti-quark to form at first a small colourless pre-hadron which rapidly 
expands to hadronic size. The time tp signifies the production time of the pre-hadron and th its later hadronisation time. 
It is argued in the text that in general tp is much less than th, and the existence of two such time scales is critical to the 
parton-hadron dynamics. 
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FIG. 2: Total and Elastic cross-sections: Comparison of data compilations from a variety of sources, CERN-HERA, UA1,4,5 
and CDF, with the fits used in LUCIFER. 
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FIG. 3: Multiplicity Distributions: Data for s^^^ = 26 GeV is compared to a histogram using the KNO scaling outlined on this 
figure. The histogram is obtained from a reasonable number of pp events, and of course is itself normalized to a probability 
distribution. Clearly the same KNO distribution is used at all energies with x = n/n. 
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FIG. 4: Multiplicity Distributions; Data for s^^^ — 546 GeV is compared to a histogram using the KNO scaling outlined on 
this figure. The strict KNO scaling seems to fail in the high multiplicity tail of the distribution shown here, but in a fashion 
which very little influences our A+A simulations. Clearly the same KNO distribution is used at all energies. 




FIG. 5: pp Pseudo-rapidity spectra; Comparison of UAl minimum bias 200 GeV NSD data 43] with an appropriate LUCIFER 
simulation. The latter is properly constrained by experiment and is an input to the ensuing AA collisions; thus does not 
constitute a 'set' of free parameters. 
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FIG. 6: A similar transverse momentum tt" spectrum from PHENIX pp [isl ] vs simulation. 




FIG. 7: Central PHENIX vr 200 GeV for Au+Au vs simulation. Curves for different choices of the production time Tp differ 
very little, since in effect the cascade effectively begins somewhat later, near 0.25 — 0.35 fm/c and continues much longer to 
tens of fm/c. Centrality for PIfENIX is here 0% — 10%, roughly for impact parameters 6 < 4.25 fm. in the simulation. 
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FIG. 8: The vr transverse momenta yields for phase I, no final cascade, vs those for the full phase I+II calculation. Clearly 
there is considerable suppression in the final cascade. Recalling that the experimentalists quote a 'direct' suppression of ~ 4-5 
for the ratio in Eqn.(l) at the highest p±, there is at the end of I an enhancement ~3, i. e. still a Cronin effect in this first 
phase. 
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FIG. 9: Pseudo- Rapidity and rapidity spectra for charged mesons and pre-hadrons at various phases of the collision simulation. 
The Gaussian fit (dashed line) to the charged pion rapidity distribution approximates preliminary BRAHMS dN'^'^/dy 
results, certainly in its FWHM, and hence demonstrates the validity of the simulation (solid line) for dN'^^ /dy. A calculation 
for dN'^'^ /dr) (dotted line) is also presented. For the charged pseudo-rapidity distributions, successive retreats, first to phase I 
of the simulation (solid circles) and then to only pre-hadrons in phase I (solid rectangles), i. e. no decays, indicates both the 
reduction in cascading participators and in the fraction of available at the earliest moments of the cascade. 
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FIG. 10: RauAu/nn- Ratio of tt" production in Au + Au collisions to that in NN at s^/^= 200 GeV for 0%- 10% centrality. The 
binary collision number selected at this centrality is Ncoii=950, while the divisor for both PHENIX and LUCIFER curves shown 
here are ones fit to the PHENIX pp data. The resultant ratio is especially sensitive to the PHENIX NN data. LUCIFER curves 
are shown for two choices A and B differing not in the Au + Au calculations but only in the fitted PHENIX NN extrapolation to 
lowest p^, where of course there is no data. The extra points for the PHENIX AuAu ratio at the four lowest p±^ GeV/c are from 
more recent analysis ^51]. Thus the general shape and magnitude of the simulation conforms well with this most recent data 
analysis, both exhibiting a rise toward smaller px. This feature is also present in the very similar 0% — 5% measurement [5ll |. 
The calculated suppression is somewhat excessive at low px • 
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